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Abstract. We compute the amplitudes for the insertion of various operators in a quark 2-point 
function at one loop in the RI' symmetric momentum scheme, RI'/SMOM. Specifically we focus 
on the moments n = 2 and 3 of the flavour non-singlet twist-2 operators used in deep inelastic 
scattering as these are required for lattice computations. 
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1 Introduction. 



Lattice gauge theories are the main theoretical tool for exploring the non-perturbative regime 
of the strong nuclear force by simulating the underlying quantum field theory which is the non- 
abelian gauge theory of quarks and gluons, (QCD). The numerical techniques allows one to 
explore the infrared region where perturbation theory becomes impractical because the value of 
the parameter controlling the expansion, which is the coupling constant, becomes large. Briefly, 
one uses the path integral formalism but with the space-time continuum replaced by a discrete 
Euclidean grid. Then one can construct Green's functions numerically, representing, for exam- 
ple, bound state particle spectra, and make measurements of the masses. Whilst overlooking 
many of the technical aspects of this procedure at this point, it is a remarkable achievement 
that the formalism is in more than solid agreement with nature. Aside from determining par- 
ticle spectra, one of the current main problems is to measure matrix elements of operators in 
the non-perturbative region. These are important in, for instance, understanding the structure 
of nucleons if one focuses on the operators relating to deep inelastic scattering. These were 
introduced originally in pQ which subsequently produced an intense industry to determine the 
operator anomalous dimensions for arbitrary moment to eventually three loops in the MS renor- 
malization scheme, [21 El [H EJ EJ [7] . Indeed there has been a large degree of progress in the 
area of measuring matrix elements of quark bilinear currents and operators and the associated 
renormalization constants by the QCDSF collaboration, [U El [TUl [III 021 [TU [14] , and others 
[13 HQ [HI H9l [SOI [SB H3l EU- However, in order to make reliable measurements and 
hence accurate predictions one has to overcome various theoretical difficulties. Aside from those 
relating specifically to the lattice, there is the problem of ensuring the results match on to what 
would be expected at high energies. In other words the Green's function depends on some refer- 
ence momentum value and the numerical simulations, in principle, make measurements not only 
at low but also at high energies. In the latter case perturbation theory is actually valid there and 
hence is a reliable complementary tool. Therefore, if the same Green's function is computed to 
several loop orders then it ought to be the case that the numerical measurements will overlap at 
large energies. Given this, it is sometimes the situation when there are accurate large loop order 
results that the continuum estimate is used to assist with normalizing the lattice measurements. 

This brief overview clouds some of the more technical aspects of the overall procedure. For 
instance, all the operators of interest undergo renormalization. Whilst this is not a major 
problem, since the formalism to carry out a renormalization has been established for many 
years now, there is the problem of the relation of a continuum renormalization to what is 
performed in practice on the lattice. For instance, the standard practice in high energy problems 
is to dimensionally regularize QCD and then to subtract the resulting divergences, which are 
manifested as poles in the regularizing parameter e, in a minimal or modified minimal way. The 
latter scheme, MS, is the main procedure primarily as convergence is improved by removing a 
specific finite part in addition to the basic poles. The main advantage of this mass independent 
renormalization scheme is that one can compute to very high loop order when the quarks are 
massless and even to a reasonable order in some cases when there are massive quarks. Whilst this 
provides accurate results for the lattice to match to, there is a technical problem to be overcome 
which is that the lattice computations are invariably in a non-minimal renormalization scheme. 
So to make a proper comparison for matching at high energy one has to convert the results to 
the same scheme. One of the more widely used lattice schemes is the RI' scheme, |25}I26| . which 
denotes the modified regularization invariant scheme. It is a modification of the RI scheme, 
[251 126] . where the essential difference is in the way the quark wave function renormalization 
constant is defined. Briefly, the difference between the RI' and RI schemes is in a differentiation 
of the quark 2-point function with respect to the momentum. As the derivative has a financial 
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cost for the lattice, the RI' scheme is more efficient and hence is the default scheme in this 
respect, |251 126j . Both these lattice schemes are mass dependent renormalization schemes and 
are a hybrid of MS and MOM schemes. By this we mean that in the main 2-point functions 
are renormalized according to a MOM type subtraction whilst three and higher point functions 
are renormalized using MS. This ensures, for example, that the RI and RI' scheme coupling 
constants are the same as that of the MS scheme, [251 12S]. Whilst introduced in [25J [26], the 
renormalization of QCD in the continuum has been studied in the Landau gauge in |27) and 
later in a general linear covariant gauge in |28j . 




Consequently, with interest in measuring matrix elements on the lattice relating to nucleon 
structure there has been a need to carry out the continuum RI' renormalization of the same 
Green's functions. These involve the low moment flavour non-singlet twist-2 Wilson operators 
which arise in the deep inelastic scattering formalism. Indeed three loop results are available in 
RI' in [28} [29| [30] . As these renormalization constants are determined by inserting the operator 
into a quark 2-point function with massless quarks then one could apply standard algorithms 
such as Mincer, [31[ I32j. to achieve high loop orders. However, from a technical point of view 
this renormalization is carried out at a point of exceptional momentum since the operator is 
inserted at zero momentum. In this configuration one is effectively dealing with a reduced 2- 
point function. It has recently been pointed out, [33], that for matrix elements relating to quark 
masses then one could have infrared issues in extracting reliable results in the non-perturbative 
region as a consequence of this momentum configuration. To circumvent this technicality an 
alternative scheme has been developed which is called the RI'/SMOM scheme, [33]. It differs 
from the RI' scheme in that 3-point functions are not subtracted at an exceptional point but 
instead at a symmetric point. This means that none of the external momenta are nullified, so 
that the potential infrared singularity embedded within the logarithms of the Green's function 
are bypassed, [33]. Whilst the one loop computations were performed in [33J, this has recently 
been extended to two loops in |34} [35] for the scalar and tensor currents. However, given that 
there is recent interest in measuring twist-2 flavour non-singlet operator matrix elements on 
the lattice for low moments, [23J, it is the purpose of this article to present the first one loop 
computations for the moments n = 2 and 3. Although lattice computations focus on the Landau 
gauge, partly as that gauge is less complicated to fix numerically, the Green's function with 
the operator insertion is gauge dependent. So in developing another renormalization scheme 
to overcome one technical problem there are potential numerical errors in measurements from 
gauge fixing issues. Instead all our results will be in a general linear covariant gauge. Whilst 
the renormalization of these Wilson operators has already been carried out for the RI' scheme, 
[281 [29| 130] . it may appear to be rather simple to follow that procedure for the latest scheme. 
This is not the case for an elementary reason. This is to do with the fact that the Wilson 
operators mix under renormalization, [36J. It is widely accepted that the flavour singlet Wilson 
operators mix among themselves and the flavour non-singlet ones do not, [H [2], [3] . Indeed in the 
original context of pQ El [3] this is the situation. However, that is only the case for the latter set 
if the Green's function containing the operator does not have a momentum flowing out of the 
operator itself. If there is a net momenta flow through the operator then they mix into a set of 
total derivative operators. Given that the RI'/SMOM scheme is at non-exceptional momenta 
then this mixing cannot be ignored. It has been studied to three loops in a practical situation 
in [36] where a similar problem for the lattice was examined but in a context which involves 
a Green's function which is gauge independent but contains two operators. Indeed the mixing 
matrix for the two Wilson operators we consider here was calculated to three loops in the MS 
scheme. Without including the mixing in [36j the operator correlation function did not correctly 
satisfy the renormalization group equation at two and three loops. 

The article is organised as follows. We review the problem and the necessary background in 
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the following section. All our results are collected in sections three and four where we give all 
the amplitudes as a function of the gauge parameter of an arbitrary linear covariant gauge for 
both the MS and RI'/SMOM schemes. We conclude in section five and include an appendix. It 
contains the bases of tensors into which all the Green's functions are decomposed together with 
the coefficients of the projection tensors which project out each specific amplitude. 



2 Preliminaries. 

In this section we recall the background to the problem and the calculational set-up. First, the 
various operators we will be considering, which are gauge invariant, are 
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where the first three operators are included for checking purposes and all derivatives, both 
ordinary and covariant, act to the right. In (|2. 1 j) S means total symmetrization in the free 
Lorentz indices and we use the same labelling and notation as [36] for ease of reference. For 
instance, at certain points in this respect we will refer to the level Wi or W3. By this we will 
mean either the specific operator with that label or the set of operators within that level which 
are additionally either dW 2 or dW 3 and ddW 3 respectively. Like [36J, it will be clear from the 
context which is meant. As discussed already one must include these additional total derivative 
operators within each level since there is mixing between the operators. Such mixings must 
be included when there is a momentum flowing through the operator insertion in the Green's 
function irrespective of the number of such included operators. For the earlier RI' scheme 
computations of the anomalous dimensions the mixing issue was not relevant since the operator 
was inserted at zero momentum, [28} l29j 130] . Given that we are considering operators with 
free Lorentz indices, whether they relate to deep inelastic scattering or not, we cannot follow 
the earlier prescription of [H EJ [3] . There the free Lorentz indices of the matrix elements were 
contracted with a null vector, A^, which projected out that part of the matrix element containing 
the divergence. The reason that we have to take a different approach resides in the fact that 
on the lattice measurements are made for various individual components of the free indices. 
Therefore, we have to take a more general approach and decompose our Green's functions into 
a basis of Lorentz tensors which have the same symmetry structure as the operator which is 
inserted into the Green's function. For the tensor current this means that the basis has to 
be antisymmetric in the two free indices since a^ v = \ [7^ , 7"] . In the case of the two Wilson 
operators each operator in the respective levels are totally symmetric in the indices and are 
traceless in d-dimensions, [lj. To be specific we have 
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where the basic operators are 



O^ 2 = ^7m^» 
= H^D^ 



(2.4) 



and Z}^ is the usual covariant derivative. These are the same definitions as used in earlier 
computations, [281 1291 130] . The total derivative operators in the same respective levels satisfy 
these same template definitions. (We have suppressed the free flavour indices but note that 
these are flavour non-singlet operators.) 

Next we recall the key points concerning the mixing of the operators in levels W2 and W3. 
First, as we are working with massless quarks there are no lower dimensional operators to be 
included and there is no mixing between levels. Next the particular choice of operators, (12.ip . 
means that whilst there is mixing the mixing matrix of renormalization constants is upper 
triangular and given by, |36| . 
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Again we avoid a clumsy index on the matrix elements by using a numerical map to the respective 
sets {WjjjoWij} and {W3, dW%, ddWz} respectively. Here the superscript indicates the level. 
Once these have been determined in a specific renormalization scheme then the anomalous 
dimension matrix is deduced from 
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where a is the gauge parameter of the canonical linear covariant gauge and a = corresponds 
to the Landau gauge. Although the operators we are considering are gauge invariant, in a mass 
dependent renormalization scheme, such as RI' or RI'/SMOM, the anomalous dimensions can 
depend on the gauge. This is why we have included the second term on the right side of (|2.8p . 
However, for the one loop computation here the leading term is scheme independent so that there 
is no gauge dependence at this order. In order to compare with the structure of our RI'/SMOM 
results later, we recall the three loop MS scheme anomalous dimension mixing matrices, [36], 
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where Q(z) is the Riemann zeta function, a = g 2 / (167T 2 ) and Nf is the number of massless quarks. 
The group Casimirs are defined by 
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where T a are the Lie group generators with structure functions f abc and 1 < a < N A where N A 
is the dimension of the adjoint representation. 

We turn now to the set-up for the particular Green's function we are interested in which is 
(^ip(p)O t fM1 ^ n (— p — q)'4 ) {q)^ and is illustrated in Figure 1. The independent external momenta 
we use are p and q and are the momenta flowing into the external quark legs. Thus there is a 
momentum of (p + q) flowing out through the operator insertion whose location is indicated by 
the circle containing a cross. In order to determine the renormalization constants for the basic 
operators in the MS and RI' schemes one chooses q = — p. However, for the RI'/SMOM scheme 
the two momenta are left unconstrained. Instead to define the symmetric point of subtraction 
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Figure 1: Graphical illustration of the Green's function, 
renormalize operators in the RI'/SMOM scheme. 
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where \i is the renormalization scale introduced to ensure the coupling constant is dimensionless 
in d-dimensions. Given this the Green's function is decomposed into a basis of independent 
Lorentz tensors, ^ n (p, q), with associated amplitude, ££,l(p, q), which is the value we 

will compute at the symmetric subtraction point, 
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Here the bracketed subscript k labels the tensor of the basis and the superscript i is the operator 
level label, (|2.ip . The explicit tensors for each level are given in Appendix A together with the 
method which allows one to compute the amplitude itself via a projection onto the Green's 
function with free indices. The same tensor basis and projection is used for each level. The 
total number of projectors, rij, is different for each level and recorded in Table 1. It is worth 
noting that the basis of projection tensors which we use for each level is not unique. They are 
constructed from the basic momentum vectors, 7-matrices and metric tensors available, in such a 
way that each final tensor has the same symmetry structure as its associated operator insertion, 
as well as being traceless in d-dimensions. However, to construct all the tensors in each set 
of projectors we have isolated all the basic one loop tensor integrals within each computation. 
These are then decomposed into their own tensor basis with their associated scalar integrals. 
For instance, 
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are the basic scalar integrals. This together with the other relevant tensor integrals are substi- 
tuted back into the main computation and the index contractions are performed. This ensures 
that no basis tensor is overlooked in the procedure of constructing the basis. Once the basis is 
established we determine the explicit projection of Appendix A. This is used to check the final 
results. Moreover, having found the tensor basis at one loop it can then be applied directly in 
a two loop calculation. 
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Table 1. Number of projectors for each operator insertion. 



In addition, for the 7-algebra we use the generalized basis of 7-matrices which has been 
introduced in [37], [38], [39] with more explicit details of their properties given in [10] 141]. Briefly, 
we define these new matrices, T^j"' 4 ", by 

r Mi...Mn = 7 [ w _ _ _ 7 Mn] (2.17) 

which is totally antisymmetric in the Lorentz indices for n > 1 and the square bracket notation 
includes the overall factor of l/n\. So, for instance, Vm\ is the unit matrix, FVq) = I, and 
a iiv _ r^. Though we will invariably retain 7^ for n = 1. One advantage of this basis is that 

tr (rg^r^p) oc 5 mn I^-^-^ (2.18) 

where JMi-Mm^i-^n j s t h e un i t ma t r ix in this T-space, [37] [381 S3 HQ] E]. As a result of this 
there is a partitioning of the projection matrix into sectors with different F( n ) as is evident in 
the examples in Appendix A. Two final points concerning the projectors are worth emphasising. 
First, our choice of basis tensors, T 7 ^ (p,q), only has meaning strictly at the symmetric 

subtraction point. Away from this point there will be a bigger basis set of tensors since then 
we would have p 2 7^ q 2 , p 2 7^ (p + q) 2 and q 2 7^ (p + q) 2 as is evident from the explicit forms 
given in Appendix A. Second, at the symmetric point the number of independent tensors in 
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each case is clearly larger than that of the asymmetric forward subtraction point considered 
in the RI' scheme. For instance, in the case of the scalar there are two tensors in the basis 
unlike the unique one of the RI' scheme. In this case the additional basis element is r?^ where 
we use the convention that the appearance of a momentum vector in place of a Lorentz index 
indicates the contraction of that index with that momentum. Clearly this object vanishes if 
either momentum is zero or one is proportional to the other, corresponding to the RI' scheme 
momentum configuration. As regards the tensor basis for each of the operators, the use of the 
generalized matrices T^"^ n is important in regarding the basis as complete. This is because 
they span spinor space in d-dimensions, |37[ I39j. with n being any positive integer. Therefore, 
it is natural to make use of them in dimensionally regularized computations. Though for the 
operators considered here n never exceeds 4. If there were more than two independent momenta 
then obviously a larger value of n would be necessary. As r^"^ n clearly form the basis of the 
spinor vector space of the tensor basis decomposition, the Lorentz vector space part of the overall 
basis for each operator is then made complete by building Lorentz tensors from combinations 
of elements of the set {77^, p^, q^, r^" Mn }. These have, of course, to be consistent with the 
symmetries of the Lorentz indices of the operator inserted in the Green's function. From this 
it is clear that our basis is complete for each operator, as there is no room to build additional 
tensors from the basic structures of the Lorentz and spinor vector spaces. 

Given that there is more than one amplitude for each operator insertion, we have to be careful 
in defining the renormalization constant in the RI'/SMOM scheme. For all the cases we consider 
here the ultraviolet divergence resides in a subset of the amplitudes which in fact contains at 
least one element except for the special case of the vector current. For V the RI'/SMOM scheme 
renormalization has to be treated separately. If, for the moment, we denote this representative 

basis tensor by the label then we define the renormalization constant for the operator O, 

^RI'/SMOM , , , 
Z , by the condition 
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where is the quark wave function renormalization constant in the RI' scheme which is 

given in [27J EH] • The reason why the value in the original RI' scheme is used has been discussed 
in |33|. [34} [35] . In determining the final renormalization constant Z^ /^MOM ^ ^ follow the 
procedure of |42j for automatic Feynman diagram computations. In other words we compute all 
diagrams in terms of their bare quantities which here are essentially the coupling constant and 
the gauge parameter. Then the renormalized parameters are introduced by rescaling with the 
already determined coupling constant and gauge parameter renormalization constants. Although 
the latter should be taken to be in the RI' scheme to the one loop order we are working any 
scheme effect will not show up until two loops. Whilst this is a standard procedure for introducing 
counterterms, the main issue here is that this rescaling from bare to renormalized quantities must 
also include the mixing of the operators. Therefore, in constructing our amplitudes, which are 
recorded in sections three and four, the matrices (|2.5j) and (j2.6[) have been included. In practical 
terms this means that the renormalization constants are found by first fixing those in the last 
row of each matrix. Then those in the next row are determined and repeated until the ultimate 
row is found. This is similar to the method used in [36] to deduce (|2.9p and (|2.10p . The 
basic reason for this literal bottom up approach is that the counterterms to be determined are 
intertwined due to the triangularity of the matrix and this is the systematic way to disentangle 
them. However, in defining the RI'/SMOM scheme here in the general terms indicated in (12.191) 
one has to be careful in any situation involving the vector current due to the underlying Slavnov- 
Taylor identity which affects this operator. We will discuss this caveat in more detail later as it 
arises not just in the case of the vector current itself but is embedded in each set W% and W3. 
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Having concentrated on the general quantum field theoretic formalism that is used, we now 
comment on the practicalities of the calculation. We use standard tools for this. All the algebra 
is carried out with the symbolic manipulation language Form, |43| . The three one loop Feynman 
diagrams are generated in electronic form by the Qgraf package, [H], with the output converted 
into Form input notation. This procedure appends indices and labels to all the fields. Various 
Form modules were then used to insert the basic Feynman rules for the propagators and vertices 
before those of the particular operator of interest. The rii amplitudes were then projected out 
by the theory given in Appendix A for successive Green's functions. These are first determined 
by constructing all the basic tensor integral decompositions akin to (12. 15H . Once determined 
we construct the tensor basis and repeat the computation by use of the projections given in 
Appendix A. The final part is to actually substitute the explicit values of the master one loop 
scalar integrals. This may have required integration by parts but for the most part the resulting 
integrals are one loop with only two propagators such as I\. These bubbles are simple to 
determine. The remaining integral is 



where S2 = (2-v/3/9)Cl2(27r/3) and Ck>(x) is the Clausen function which was evaluated in [45] . 
We use dimensional regularization in d = 4 — 2e dimensions. Given these ingredients we have 
been able to determine all the one loop amplitudes for the set of operators f)2. If) where we have 
repeated the calculation for S, V and T as an elementary check on our programmes. We correctly 
reproduced all the one loop expressions given in [33^ [34"1 [35] for both anomalous dimensions and 
amplitudes. 



(2.20) 



3 Quark currents. 



Having concentrated on describing the background to the problem and the methodology of the 
computations, we turn to the mundane task of recording the explicit amplitudes. In this section 
we do this successively for the scalar, vector and tensor currents. In [33[ [34"1 [35] the anomalous 
dimensions were computed at one and two loops for the scalar and tensor cases in the RI'/SMOM 
schemes. In addition the one loop amplitudes in the MS scheme were given for each of the 
three operators in [33]. In the appendix we give the explicit mapping between our amplitudes, 
T,n,\(p, q), for the quark currents and those defined in [33] in order to compare. Though we note 



that we are in full agreement with the results of [33J. Since lattice computations will measure 
various directions and then extract estimates for the overall Green's function it seems appropriate 
to provide the amplitudes for the MS scheme, as well as for RI'/SMOM, as ultimately the former 
is the reference scheme one will map to. Equally the explicit MS expressions will be useful for 
lattice groups who convert to MS first before comparing to the high energy expressions rather 
than work in the RI'/SMOM directly in order to do the matching. This seems appropriate since 
it will be noted later that there is not a definitive RI'/SMOM scheme especially when operators 
have Lorentz indices such as the Wilson operators. Moreover, whilst providing results for both 
schemes may appear to introduce a degree of redundancy, because at one loop there is overlap 
in the actual expressions for the amplitudes between the schemes, this is not preserved at two 
loops, |46] , Also for the Wilson operators the results in both schemes will be useful to ensure 
that we have preserved properties of the sets of operators in each level between schemes. For 
the scalar current we have the MS result at one loop 
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To make contact with other work in this area, [33], we note the relations 
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where ip(z) is the derivative of the logarithm of the Euler T-function and Cq{uj) was the function 
used in [33] to interpolate between the symmetric and exceptional momenta scheme choices. 
Whilst this case has already been analysed in [331 EH EE] we discuss it here as it raises several is- 
sues with regard to defining the subsequent RI'/SMOM scheme renormalization constants for the 
Wilson operators. For instance, it would appear that the definition of the scheme can be given 
in several ways. Given this form for the Green's function, one could either define the RI'/SMOM 
scheme operator renormalization constant by absorbing the finite part of the 1 direction. Alter- 
natively one could take the projection of the Green's function by some tensor, such as that of the 
Born term, and then absorb whatever the finite part emerges into the renormalization constant. 
For the scalar case this would actually give the same renormalization constant because of (|2. 18j> . 
However, for other operators there appears to be a degree of freedom in how one can actually 
define the scheme. Though, of course, one would retain the fundamental criterion of (|2.19p to 
ensure there are no corrections beyond the leading order. As these methods are equivalent here 
we note that we reproduce the scalar current results of |33j for the renormalization constant 
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where we use the notation that all expressions will be in RI'/SMOM unless explicitly indicated 
to be in MS. The resulting amplitudes are 
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where channel 2 is the same as the MS expression. However, this will not be the case at two 
loops, 



The situation for the vector current is more involved partly because of the increase in am- 
plitudes with one free Lorentz index. First, the MS amplitudes are 
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Here the renormalization constant for the current is 



rV 



MS 



1 + 0(a 2 



(3.5) 



(3.6) 



This value derives from the fact that the vector current is a physical operator and undergoes 
no renormalization to all orders in perturbation theory. Moreover this is a consequence of 
the Slavnov- Taylor identity. Turning to the RI'/SMOM situation, it is tempting to define the 
renormalization constant in this case by either the projection by the Born tensor or by taking 
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the coefficient of the channel 1 amplitude. This would lead to a renormalization constant in the 
RI'/SMOM scheme with a non-zero finite part. However, this would be inconsistent with the 
Slavnov- Taylor identity. In addition as the MS renormalization constant is unity to all orders 
then this implies, [IT], that the RI'/SMOM renormalization constant is already determined and 
equivalent to (13. 6ft . Therefore, we have for the vector current renormalization 
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(3.7) 



In order to see that this is equivalent with the Slavnov- Taylor identity one can contract the 
Green's function with the vector (p + g) M as this is the momentum flowing through the inserted 
operator. This procedure corresponds to the renormalization of the divergence of the current. 
In terms of (|3.5p the combination of amplitudes in this contraction proportional to $ is 
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A different combination determines the coefficient for the piece involving ^ but with the same 
result. The fact that there is no 0(a) correction is because we have renormalized the quark 
2-point functions in an RI' scheme and the wave function renormalization constant defined in 
this way leaves the 2-point function as unity to all orders. Hence the Slavnov- Taylor identity is 
satisfied. We have checked that this is also the situation in the MS case and, moreover, this has 
been extended to two loops in RI'/SMOM, [46], where it was verified that the Slavnov- Taylor 
identity was consistent in that case too in keeping with the general argument given in [33]. 
Crucial to this analysis is knowledge of the full basis of tensors into which the Green's function 
can be written. Although the channel 6 amplitude appears to be decoupled due to the rv n ) 
basis we use, omitting it would lead to an inconsistent operator renormalization. Finally, in 
considering the Slavnov- Taylor identity in this way it would have become clearer if a different 
combination of tensors was used in the basis rather than the ones given in Appendix A. In other 
words we could have used a basis where all but one tensor in the basis vanished when contracted 
by (p+q)u whence the presence of the Slavnov- Taylor identity would have been explicit. Finally, 
we note the RI'/SMOM scheme amplitudes after renormalization are 
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(3.9) 



As this particular Green's function is symmetric under swapping the two independent momenta, 
then two pairs of the amplitudes are equivalent. The different values for E^(p, q) in both 
schemes is due to the fact that the finite parts of the quark wave function are not the same in 
each scheme. 

For the tensor case we provide the same information as the previous cases. First, the MS 
amplitudes are 
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(3.10) 



The earlier comments concerning how one actually defines the RI'/SMOM scheme apply here. In 
the finite part of the tensor current renormalization constant was defined by contracting 



the Green's function with and ensuring that the resulting expression had no a dependence 
after renormalization. Following that with our amplitudes we find exact agreement with the 
divergent and finite parts of the same renormalization constant as, [33], 
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As in the vector case with this projection to define the renormalization constant the T^) and Tu\ 
sectors do not contribute but are part of the full Green's function at the symmetric subtraction 
point. Indeed we have checked that the full expression for the two loop amplitude, [36], leads 
to the same renormalization constant as [35]. However, if one were to follow an alternative way 
of defining the RI'/SMOM renormalization constant by merely using the coefficient of channel 
1 which is where the divergence resides then one would have another expression for Z T which is 
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This alternative definition is of course dependent on the choice of basis tensors. One feature 
of it is that the numerical value of the finite part in the Landau gauge is significantly smaller 
than the corresponding part of (|3.11|) . Indeed with a judicious choice of the projection basis it 
might be possible to have a rapidly converging series for the conversion function. However, that 
requires a higher loop analysis, [35]. Returning to the original RI'/SMOM scheme definition of 
, [35] the amplitudes are 
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(3.13) 



Similar to the vector case there is a clear symmetry here because the Green's function is sym- 
metric under the interchange of p and q which is responsible for the two pairs of equivalent 
amplitudes which is clear from the explicit tensor definitions given in Appendix A. 



4 Deep inelastic scattering operators. 



The situation for both moments of the flavour non-singlet twist-2 Wilson operators is more 
involved due to the operator mixing issue as well as the issue with the Slavnov- Taylor identity 
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noted previously. First, the MS amplitudes for n = 2 are 
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1.1) 



where the operator superscript label here corresponds to the row of the matrix with that operator 
on the diagonal. Unlike the previous two cases there is now no symmetry for the Green's function 
itself when the original operator is inserted. This is because the covariant derivative in the 
operator only acts on the quark and not the anti-quark. So swapping the external momenta in 
the Green's function is not a symmetric operation. By contrast for the associated total derivative 
operator this symmetry is still valid which is why there are equivalences between various pairs 
of amplitudes. Moreover, the actual expressions for the tensors labelled 3 and 5 are proportional 
to those respectively labelled 2 and 3 of the vector case. This is not unexpected because the 
total derivative operator associated with this moment is effectively the vector current in disguise 
at the higher level. That the expressions are not precisely the same is due to the fact that one 
has an extra Lorentz index present at this level so that the projection coefficient into the basis 
will not have a completely parallel value. Although there is no unique basis for the projection 
tensors this agreement at one loop is a check on our calculational setup particularly since with 
another choice of tensors this proportionality could have been hidden. For the renormalization 
constants the relation between the divergences has already been noted at three loops in the MS 
scheme in |36|. 



Now that we have all the finite parts at the symmetric subtraction point we can define the 
RI'/SMOM scheme renormalization constants. Similar to the vector current this requires care in 
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the case of dW 2 - This operator is the total derivative of the vector current whose renormalization 
constant is already determined in all schemes as it is a physical operator. However, only when 
one takes the contraction of the two free Lorentz indices of dW?, does the divergence of the vector 
current emerge. Therefore, to have a renormalization consistent with Z v we have to define Z^ 2 
to be unity in the RI'/SMOM scheme. In the MS case, to contrast with (13.8p . the appropriate 
combination of amplitudes for the piece proportional to $ gives 
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+ -C F aa + 0(a 2 ) (4.2) 



whence it can be recognised that the Slavnov- Taylor identity is preserved. An alternative choice 
gives the $ contribution but it has the same result as (|4.2p . This combination is deduced by 
contracting with (p + q)fi(p + q) u to ensure each term in the definition of the operator involves 
the divergence of the vector current. However, for the other operator of level W 2 there is no 
underlying Slavnov- Taylor identity as the contraction of the two free Lorentz indices is not 
the divergence of a physical current. Therefore, there is no constraint on how to define the 
remaining two elements of the W2 renormalization constant matrix in RI'/SMOM. Indeed in 
some sense there is an infinite choice. However, for our purposes here we have chosen to define 
these renormalization constants by ensuring that there are no 0(a) corrections to channels 1 and 
2 which both contain the divergences in e. The former contains the off-diagonal counterterm of 
the mixing matrix whilst the latter contains both counterterms from the first row of the matrix. 
Clearly one fixes the off-diagonal one first. Therefore, we have the W% matrix of renormalization 
constants 
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where in defining Z^ 2 the right hand side of the analogous expression to (|4.2p we have ensured 
that there are no 0(a) corrections. Consequently, the RI'/SMOM amplitudes are 
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(4.4) 



At this loop order the finite parts of a substantial number of the amplitudes are the same as their 
MS counterparts. However, we do not expect this to extend necessarily to two loops. Finally, 
we note that there is a relation between various amplitudes 
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which may still be valid at higher loop order. 
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For the next moment, n = 3, the situation is of course more substantial since there are four- 
teen different basis tensors and three operators which mix. First, we record the MS amplitudes 
to one loop are 
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Whilst there are more amplitudes some features are similar to level W 2 such as equivalences 
between certain amplitudes for ddW^ and a proportionality with various amplitudes in V and 
dW 2 . As this is a level higher there is also a proportionality of certain dW^ amplitudes with 
W 2 . For instance, the channels 2, 3 and 4 respectively of V, dW 2 and ddW$ are proportional 
as well as the respective channels 3, 5 and 7. The same is true for W 2 and dW% for the two 
respective pairs of channel 3 and 4 as well as 5 and 7. This is of course not unexpected since the 
operators are successive total derivatives of the lower level one. However, the renormalization of 
ddWs in the RI'/SMOM scheme is already predetermined for the same reasons as dW 2 . Though 
in the case of ddW% one has to contract with (p + q)^(p + q) v {p + q)a to produce the relevant 
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divergence of the vector current. Using (p + q)^rj ua as an alternative nullifies the operator as 
it is constucted to be traceless and gives nothing non-trivial. With the former contraction this 
produces the relation between the amplitudes involving $ is 
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An analogous relation produces the piece involving $ which also has no 0(a) piece consistent with 
the RI' choice for the quark wave function renormalization. Clearly this combination ensures 
consistency with the Slavnov- Taylor identities. 

Again the choice of how to determine the remaining renormalization constants within the 
RI'/SMOM scheme ethos is relatively free. The point of view we take to do this is to build on 
the previous level. Though there is no reason why one should necessarily make this the definitive 
choice. By building on the previous level we parallel the way we ensured that the hidden vector 
current renormalization constant was determined. As dW% is the total derivative of W% we 
choose that level of the renormalization constant mixing matrix to be the same renormalization 
constants as the first row of the W2 matrix. This leaves the first row of W3 and we then define 
this by ensuring that after renormalization there are no 0(a) corrections in the tensor basis 
which originally had a divergence in e. As there is mixing within the counterterms we have to 
determine these in sequence which was the same as the MS scheme, [36]. Therefore we have 
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Once these have been determined the RI'/SMOM scheme amplitudes are 

^ E$(p,g) = 0(a 2 ) , 
+ 0(a 2 ) , 
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As with the lower moment there are similar relations to one loop. First, there is an obvious 
parallel with the vector and total derivative operator for n = 2 with the double total derivative 
of the third moment which was noted in [36]. In addition there are cross connections with 
the central row of the mixing matrix. Specifically, the amplitudes labelled 6 and 8 of W2 are 
proportional to 8 and 11 of dW% respectively for the same reasons as before. Further, there is 
another apparent relation within the level since 

Eg^-fog) - 2E*™ r »(p,g) + 47 3 (P,<?) = 0(a 2 ) (4.10) 

in addition to the reflection of the one already noted for moment n = 2. As an aside we note 
that in highlighting relations between amplitudes, we have concentrated on those which are true 
for all values of a. There appears to be other relations which are only valid to one loop in the 
Landau gauge and since the case for regarding these as significant is diminished because of the 
apparent gauge dependence we do not draw attention to them. 

As our computations have been at one loop the conversion functions that are used to convert 
from the MS scheme to the RI'/SMOM scheme are in effect the finite parts of the renormalization 
constants themselves. See, for example, [IT]- Therefore, in order to appreciate the magnitude of 
the one loop contribution it is a straightforward exercise to evaluate the finite parts numerically. 
It turns out that the correction increases in value with the operator moment, in parallel, for 
example, with the numerical value of the one loop anomalous dimension coefficient. If we take 
the Landau gauge case and compare the finite parts, the moments n = 2 and 3 are respectively, 
3.8436 and 6.1839 where we use the (11) element for the latter two or equivalently the diagonal 
of the W3 matrix. This is with ignoring the common colour factor, Cf, and coupling constant, a. 
Therefore, for higher moments it would seem that one would require higher loop corrections in 
order to have a more reliable convergence of the perturbative series for the conversion functions. 
Comparing with the RI' scheme results of [281 [29j [30] the respective numbers are 3.4444 and 
5.9444. Unlike the situation with the mass operator the RI'/SMOM scheme corrections for W2 
and W3 are slightly larger. Though in some sense this is not a fair comparison because of the 
mixing. The RI' scheme is based on a specific momentum routing in the Green's function which 
is blind to the off-diagonal matrix element of the renormalization constant matrix. Equally the 
convergence of the RI'/SMOM series could be improved by choosing a different combination of 
the finite parts to define the actual renormalization constant since there appears to be a large 
degree of freedom in doing this. Only a higher order computation could give insight into this, 



5 Discussion. 



We have determined the renormalization constants at one loop for the twist-2 flavour non- 
singlet operators with moments n = 2 and 3 in the RI'/SMOM scheme and have taken account 
of operator mixing. In the Landau gauge the correction increases in value with the operator 
moment. So it would appear that the series for the conversion function will converge slower 
at high moment. Of course, a one loop computation is not sufficient to make a definitive 
statement since the higher loop corrections may produce an improvement. This is currently 
under investigation, |46] , Equally it would be useful to see how the lattice measurements in this 
symmetric scheme for the deep inelastic scattering operators compares with the same analysis 
in the earlier RI' scheme. Though it may be the case that not all the various tensor projections 
would give a clear accurate numerical signal. Finally, whilst the RI'/SMOM scheme is designed 
in order to circumvent infrared problems, one is still restricted to the Landau gauge due to 
the fact that the Green's function of Figure 1 is a gauge dependent quantity. Whilst this is 
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not a problem for the high energy regime, it may be the case that in the infrared there are 
additional complications with Gribov copies and hence the definitive measurement of operator 
matrix elements. 



Acknowledgement. The author thanks Dr. P.E.L. Rakow and Prof. C.T.C. Sachrajda for 
useful discussions and especially the former for a careful reading of the manuscript. 



A Projectors. 

In this appendix we record in succession the basis of projection tensors we have used in each 
of the various cases. These were denoted earlier by the general notation Vfy^ ^ n (p,q). The 

matrix, AIL, in the second part of each subsection contains the coefficients associated with 
each basis projection tensor, which allows one to project out the amplitudes H9^(p, q). It is 
constructed by first determining the matrix 



A? 



A7 



(k) n\...Hn 



l |Ul...jU„ 

(I) 



(p,q) 



(A.l) 



where there is no summation over the label i and k and I index the projection tensors. The 
matrix J\fh is symmetric in k and I. Given that the elements are contractions of Lorentz tensors 
in d-dimensions then the elements of A/Jy are polynomials in d. Finally, AIL is the inverse of 
M^i an d then 



Sg(p,g) = Mi^-^&q) ((mo; i ...» n A-p-qMq) 
where again there is no summation over the level label i. 



(A.2) 



A.l Scalar. 



1 



P(i)(P,9) = r (0) , V {2) {p,q) - ^ (2) 
12 ^ -4 J 

In order to compare with [33], the mapping between the amplitudes is 

Ef 1} (p,5) = - A s , Ef 2) (p,g) = - 2C S 
where the notation of the right hand side is that of 



(A.3) 
(A.4) 

(A.5) 



A.2 Vector. 



7/i . "P(2)n(P>l) 



2 ( 3 ) 



(A.6) 
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For comparison with |33j the amplitudes are related by 



s Ji)(p> q) 
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(A.7) 



A v + 2B V + -C v + -By , E ( v 2) (p, g ) = E ( v 5) (p,g) = 2By , 



2 



) S (6)(^'^) = Cy - Dy . 



(A.8) 



A. 3 Tensor. 
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(A.9) 



To save space for the projection matrix in this case, we have partitioned the 8x8 matrix into 
four submatrices. We have 
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In the tensor case the relations with the amplitudes of |33j are 



s fi)(P>tf) 
£(3)(P>9) 
s fs)(P^) 
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(A.10) 
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(A.12) 
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A.4 Wilson 2. 
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(A.13) 



Similar to the previous case we have subdivided the 10 x 10 matrix here into four submatrices 
in order to ease the presentation. We have 
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A. 5 Wilson 3. 



^2 Vl^PvPa + luPaPfi + laPfiPu] + ^ 2 ] ^ ^ ua + ^VvfJ, + 
- ^ + 2 ^ 2 Wf"^ + + ^^"J > 

~2 [7/^<?<r + IvPaQti + IcP^v + l^vPa + 7i>9<tP/x + 7(7?^] 

A 1 



23 



2i> 



[d _|_ 2]/i 2 W\wVo i VvaPfi + VofiPv] i 



[d + 2] 



2^ 



[d + 2]^ 



J^P/xPuPa + ^ ~ 2 j^ 2 WvPa + TWf>M + VafiPu\ 



[P^PuQa + PnQvPa + q^PvPa] 



[d + 2]^ 



i> 

— \pp,qvq<j + q^Puqa + q^qvPa] 
i> 

i i> r i 

-j^q^qvqa + ^ - 2 ^ 2 h MI /9<T + i)voqix + , 

■J^PfiPiyPa + ^ ^]^2 + VvaPii + Va^Pu] , 

"4 tW^a + Pf^q^Pcr + q^PuPa] 

A* 



[d + 2]/x 2 
\ 

4 \p ^q^qo' + q^Pvqa + q^Pa] 



[d + % 2 



+ 
1 

M 3 



r (3) m PuPa + T( 3 ) upq Pfj,Pa + T( 3 ) apq p^p u 
1 



[d + % 2 



^(3) npqVvcr + ^(3) vp qVficr + T( 3 ) apq ri^ v 



T(3) npqPuqa + T( 3 ) upqP^cr + I\ 3 ) apq p^q u 
+ r (3) m quPa + T( 3 ) upq q^Pa + T( 3 ) apq q^p u 



r '(3) fipqVua + ^(SjvpqTlno + ^(3) apq^iiv 



[d + 2]/i 2 



^(3) npqQuQa- + T( 3 ) vpqq^qa + T( 3 ) apq q^q v 



F(3)HpqVvcr + ^ (3) vpqV^a + 1^(3) apqVp.f 



[d + 2]^ 



(A.15) 



24 



Here the partition of the 14 x 14 matrix is into ten non-zero submatrices. We have 
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